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Abstract 

Higher order renormalons beyond the chain of one-loop bubbles are discussed. A 
perturbation method for the infrared renormalon residue is found. The large order 
behavior of the current-current correlation function due to the first infrared renormalon 
is determined in both QED and QCD to the first three orders. 



1 Introduction 



Perturbation theory in field theories is generally plagued by the rapidly growing coefficients, 
which cause the series in weak coupling to be asymptotic. Classical solutions, instantons, 
cause the perturbative coefficients to grow as n! for large n, where n is the order of perturba- 
tion, and so do certain subsets of Feynman diagrams, renormalonsQ. Some of the properties 
of the renormalons are discussed here. 

A chain of the one-loop bubble diagrams in a photon propagator in massless QED is an 
example of the renormalons (Fig.l). An exchange of the one-loop Gell-Mann-Low (GL) effec- 
tive charge gives a contribution of n\ for an ultraviolet renormalon, in which the momentum 
flowing in the propagator is large compared to the renormalization scale (k 2 / /x 2 ~ exp n), and 
(— l) n n\ for an infrared (IR) renormalon, in which a soft momentum (k 2 /fi 2 ~ exp (— n/2)) 
flows in the propagator. 

The actual form of the large order behavior due to an infrared renormalon is generally 
given by 

Kn\n u b^ n (l + 0(l/n)), (1) 

where v and bo are renormalon-specific, known constants. The coefficient K is an all-order 
quantity It depends not only on the one-loop renormalon mentioned above but also on 
an infinite set of higher order renormalons, and so determining it is nontrivial. 

However, it should be emphasized that K is calculable, at least perturbatively. For 
example, if we had calculated the series to a very high order, then eq. ([I]) implies that we 
could extract the coefficient to an accuracy of 0(l/n). Therefore there must be a convergent 
sequence for K, with K being its limit, associated with the perturbation of the amplitude 
in consideration. The main purpose of this paper is to present such a sequence for the first 
IR renormalon in the Borel plane. 

The precise calculation of the large order behavior is important, besides its theoretical 
interest, because it could play an essential role in an effort to reconstruct the true amplitudes 
from the perturbation theory. The large order behavior due to the IR renormalons in non- 
abelian gauge theory arises from the imaginary part of the nonperturbative effects, vacuum 
condensations, and so a precise calculation of the large order behavior gives detailed infor- 
mation on the imaginary part of the nonperturbative amplitudes, which could be essential 
in understanding the full amplitude. For a recent consideration in this direction one may 
refer to [0]. 

This paper is organized as follows. In sec. 2, we discuss in QED the higher order 
renormalons beyond the chain of one-loop bubbles, and show in detail how the large order 
behavior gets contribution from the higher order renormalons. In systematic 
method of summing those higher order renormalons is discussed, and the renormalon residue 
to the first three orders is given. In sec. 6, we discuss the calculation of the large order 



Renormalons also denote the singularities in the Borel plane. 
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behavior in QCD using the analytic property in Borel plane, and give the large order behavior 
to the first three orders. In sec. 7, the scheme dependence of the large order behavior is 
discussed. 




Figure 1: One-loop renormalon. 



2 Higher order renormalons 

We first review how a chain of the one-loop bubble diagrams gives rise to factorial growing 
coefficients, and then show that the large order behavior of perturbation theory is an all 
order property by giving an estimate of the higher order renormalons. 

The Green's function we consider is the electromagnetic current correlation function in 
QED in the Euclidean regime 

r t o n(a(u),^) 

i j e %qx < Mx)j„(0) > d 4 x = (q„q„ - q 2 g, v ) ^ Q \ (2) 

where Q 2 = —q 2 > 0, and 

f{x)=^{x). (3) 
For large order behavior, it is more convenient to consider the renormalization scheme and 
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scale invariant quantity D defined by 



D(a(fx),fx 2 /Q 2 ) = g s 







<9Q 2 



n(a(//),&) 



-Q 2 



d 



dQ* 



a(/i)=0 



E a n (fi 2 /Q 2 )[a(fi)] 



n+l 



(4) 



n=0 



A single exchange of the GL effective charge gives rise to the first IR renormalon sin- 
gularity in the Borel plane. It is also generally assumed that the leading residue of the 
first IR singularity can be completely determined by a single exchange of the GL effective 
charge, which implies equivalently that the large order behavior due to the IR renormalon 
can be determined by a single exchange of the GL effective charge ( Fig.l ). For the other IR 
renormalon in the Borel plane, it is similarly believed that their residues can be saturated 
by the multiple exchanges of the GL effective charges. The analysis for IR renormalon in 
nonabelian gauge theory using operator product expansion supports this assumption [§, |J. 
With this assumption for the large order behavior for the first IR renormalon, D may be 
written as 

D(a(v)^ 2 /Q 2 )= I f(k 2 )*(k 2 )dk 2 (5) 
Jo 

where &(k 2 ) denotes the GL effective charge and 



f(k 2 



8vr 3 g 4 



for k" 1 



(6) 



with Cf denoting the charge of the fermion ip. This infrared limit of f{k 2 ) can be easily read 
off from the coefficient of F 2 y term in the operator product expansion of the current product 
ineq. (g). 

To see the n\ growth of the perturbative coefficient from the chain of the one-loop bubbles, 
we may substitute a(k 2 ) in eq. (|5]) with its one-loop form 



a(k 2 



at(fjL) 



1-Mfi) ln(^) 



E 

o 



a(/x) 



n+l 



(7) 
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where (3$ is the first coefficient of the (3 function, to obtain 



a 



•n 




( Po\ 
16tt 3 Q 4 \ 2 J 



n\(l + 0(l/n)), 



(8) 



where 




(9) 



For large n the leading contribution to the integral comes from the kinetic region k 2 ~ 
\j? exp(— n/2), and thus the leading large order behavior is independent of the upper bound 
of the integral. 

Let us now consider the effect of higher order renormalons on the large order behavior of 
perturbation. First we introduce some definitions. In the following the vacuum polarization 
diagrams are assumed to include two external photon propagators. An irreducible renor- 
malon is defined by replacing all photon propagators in an irreducible vacuum polarization 
diagram with chains of the one-loop bubbles. Similarly, reducible renormalons are defined 
by replacing all photon lines in reducible vacuum polarization diagrams with chains of the 
one-loop bubbles. Thus for every vacuum polarization diagram there are corresponding 
renormalons. 

We assign an order p, and the number of reduced photon propagators, q, to each irre- 
ducible vacuum polarization diagram by 



where n A is the number of photon propagators and denote the number of irreducible 

vacuum polarization subdiagrams and the number of the one-loop bubbles respectively. A 
reduced photon propagator is simply a chain of an unspecified number of one-loop bubbles. 
The same p and q of an irreducible vacuum polarization diagram are defined as the order 
and the number of the reduced photon lines of the corresponding renormalon. For example, 
the order and the number of reduced photon lines of the one loop renormalon in Fig.l is 
p = 0, q = 1. Other higher order renormalons may be similarly characterized by the p and q. 
For reducible renormalons, the highest order of the irreducible subrenormalons of a reducible 
renormalon is defined as the order of the reducible renormalon. Some examples of the higher 
order renormalons are given in Fig. 2. 



p = n A -n L 



- 1 



Q = 



(10) 
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P=3 



Figure 2: Examples of higher order renormalons. Dashed lines denote chains 
of one-loop bubbles, and (p, q) denote the order and number of the 
reduced photon propagators respectively. 



Before we discuss the effect on large order behavior of higher order renormalons in general, 
let us take some specific examples of low order renormalons, and see their contribution to the 
large order behavior. This exercise is very instructive and gives an insight to more complex 
renormalons. 

It was first noticed by Grunberg JIJ, and diagrammatically by Mueller ||, that the large 
order behavior of perturbation is an all order property ( see also ||). The following argument 
is motivated by Mueller's observation. Let us consider the order one renormalon in Fig. 3. 
The coefficient a n (t) of a(/j 1 ) n+l due to this renormalon in the perturbation of the GL effective 
charge in a(fi) is given by 



^(t) = £(-iii(f)) B - 



n 



(r-l), 



(t) (n — r) 



(11) 



r=l 
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Figure 3: Order one renormalon. Dashed lines denote chains of the one- loop 
bubbles. 

where IT r is the r-loop vacuum polarization function 

CO 

U{a^),t) = y £U r {t)a(n) r . (12) 

r=l 

The powers of 111 in eq. flTT| ) obviously come from the one loop bubbles in the external 
reduced photon propagators, and the factor (n— r) accounts for the (n—r) possible location of 
the fermion loop with the internal reduced photon line. 11^ denotes the terms proportional 
to in n r+1 . 

The general form of can be deduced by considering the following renormalization 
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group equation for IT(a(/i),t) 



2 d 

^ 1-2 



dp? 



1 1 



(13) 



which comes from the renormalization scale invariance of the GL effective charge a(fc 2 ) 



a(a(n),t) 



at(fj) 



l + n(a(/i),t)" 



(14) 



Putting the perturbative form of II in eq. (PD into eq. (0) we get the recursion equation 



d n 
t— Il„ + i = -f3 n + J2(m- l)/3 n _ m II m , 

^ m=2 



for n — 2, 3, 



(15) 



with 



ITi = -Ajlnt + p! 
n 2 = -/3ilnt + p 2 . 



(16) 



Here /5 m are the coefficients of the (3 function defined by 



/3(a(/i)) = /I 



<i/i 2 



m=0 



(17) 



and are constants. 

Solving the recursion equation we have 



(18) 



Substitution of this into eq. (|TT|) gives 



n-l 



S B (t) = X;(A)In(*)) 



ra-2 



r=l 



— ln(t) - p 2 



(n — r). 



(19) 



Here we kept only the ln(i) term in H\(t) for simplicity, and the effect of the constant term 
in H\(t) will be discussed shortly. Since a factor of (lnt) n in the integrand in eq. @ would 
give rise to 

(2°) 

a n (t) gives the following large order behavior 



e 2 ff x 4 



0o 



16vr 3 Q 4 



7?! 



-2/5 



i(lnn + 7 S -l)-2§ 



ft 2 



(l + 0(l/n)) (21) 



which is comparable to the one-loop renormalon contribution. Here je is the Euler constant. 
Going back to eq. (|ITD , expanding the factor 



-III 



.n— r— 1 



(0o lnt-pi 



,n— r— 1 



(22) 



it is easy to see that every term in the expansion proportional to 



for i <C n 



(23) 



also contributes to the large order behavior. Thus the the inclusion of the constant term in 
Hi(t) would modify eq. (^) into a form 



0o 



01 



% [{lnn)h(p 1 ) + h'{p 1 )) + ^h"{p 



ft 



;i + 0(l/n)) (24) 



where 



(25) 
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Here hi are calculable constants. Note that this series runs to an infinite order in the limit 
n — ► oo. h'(pi), h"(pi) are similarly defined in a series form. 



(DC.) • C.)(.)CD » 



Figure 4: Chains of the order one renormalons. 

Using a similar method it is now easy to estimate the large order behavior from other 
renormalons. For example, it is straightforward to check that the chains of the order one 
irreducible renormalon in Fig. 4 along with the one-loop renormalon contribution in eq. (§) 
exponentiate the ln(n) term in eq. (|2"I|) to give the n u factor in the large order behavior with 



v = ( 26 ) 
Po 



which agrees with the well-known result for the first IR renormalon [|J . 

Now to discuss the higher order renormalons in general, consider an irreducible renor- 
malon of order p and reduced photon propagator number q. The large order behavior due 
to this renormalon in the infrared regime is given by 



n-l 



n lb _L 

/ t £ (fl,lnf) n - r - 1 \n^ ii (t)]{n-r)dt. (27) 

J o 



r=q-2 



Here we picked up only the ln(i) term in Tli(t) as before. The general solution of the recursion 
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equation in eq. ( |15|) may be organized in terms of p and q as 

n r+ i = E^tE(i«r w E ?knw ( 28 ) 

qr'=l p=l fc=0 {mi,mi} «=2 

where 

= g - 2 (29) 
and mi,mi are nonnegative integers that satisfy 



Ei>2 m i + m i = q' - p + 1 

Ei> 2 K-+wij)-(*-i)=p- (so) 

Here is the constant term of Ilj(t). The Il^' ?+2 ' ) (t) of order p renormalon is then given by 

i p+i 

n^ + 7 +2) (t) = ^ +2 E(int) r - p - fc+1 E cg^nffitf 14 - (si) 

fe=0 {rrii,mi} t=2 

By solving the recursion equation in eq. fll5|) explicitly for the several low order diagrams, 
and considering the form of the large order behavior in eq. ([I]), it is not difficult to convince 
oneself that asymptotically 

C r { ^}- rP+k ' 2 9(^r), forr^oo. (32) 

Here the function g(x) is a polynomial of degree at most — 1 for m-i > 2, and causes the 
n v term in the large order behavior. 

Substituting eq. ([H]) into eq. (B7|) and ignoring the logarithmic dependence in g, we get 



a,, 



* 2 Lj2j(p,k) £ U^-\PT\ (33) 

k=0 {mi,m,i} i=2 
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where J(p, k) is a function of p and k. This shows that the irreducible renormalons of all 
order contribute to the large order behavior. Further note that the large order behavior 
due to an irreducible renormalon of order p involves only the coefficients of the vacuum 
polarization function and the (3 function to (p + l)-loop order. 

These are also true for the reducible renormalons. As we chain more irreducible sub- 
renormalons into a reducible renormalon, more powers of hit in a n (t) due to this renormalon 
are being lost, resulting in suppressed integral in t, but this suppression is exactly compen- 
sated by the larger combinatoric factor caused by the more possible locations in putting the 
the subrenormalons, giving a large order behavior comparable to those from the irreducible 
renormalons. Therefore, all reducible renormalons also contribute to the leading large order 
behavior. Also since a n (t) for a reducible renormalon with m irreducible subrenormalons is 
proportional to 

TT(ri-gi+2)- rr (r 2 -g2+2) Tr(r m -q m +2) (o A \ 
ll n+ l iJ -r2+l ' ' ' LL r m +l > \ 0q: ) 

with each Hr.+i t+2 ^ coming from the subrenormalons, it is obvious that the large order 
behavior from this reducible renormalon depends on the vacuum polarization function and 
the (3 function to (p + l)-loop order. 

The inclusion of the constant term in Hi(t) in eq. ( f27|) would modify the large order 
behavior in eq. fl33|) in a similar fashion as in the example of the order one irreducible 
renormalon in eq. Q23|)- With pi included, each term in eq. will be multiplied by a 

series in p\ in the form of eq. (|25|), with the coefficients hi now depending on each particular 
term. 

Though this analysis of diagrams is very helpful in understanding general higher order 
renormalons qualitatively, it seems difficult, or at least inconvenient, to systematically cal- 
culate the higher order renormalons using this technique. We need a more straightforward 
approach for systematic evaluation of the higher order renormalons. In the following sections 
such an approach is discussed. 



3 Borel transform of GL effective charge 

The problem of determining the renormalon residues of D(a(/j,), p 2 /Q 2 ) eventually reduces to 
finding the Borel transform of the GL effective charge. We give in this section a perturbation 
method for the Borel transform of a(fc 2 ). We are going to introduce a scheme and scale 
independent coupling, and then write the GL effective charge in terms of the coupling in a 
form that is particularly convenient for Borel transform. 



11 



The Borel transform of D (a(fi) , fi 2 / Q 2 ) is defined by 



D(a(fi),iJ 2 /Q 2 ) = f°exp( ^—)D(b)db. (35) 

Jo \ a[u) ) 



With the perturbative series of _D(a(/i), fi 2 /Q 2 ) in eq. 



00 a 



D(b) = £ (36) 



n=0 U - 



The large order behavior of the form in eq. (|1|) causes a singularity (renormalon) in D(b), 
and conversely the singularity in D(b) determines the large order behavior. Thus by study- 
ing D(b) near the renormalon singularities we can determine the large oder behavior of 
D(a^)^ 2 /Q 2 ). 

Let us first consider the renormalization group equation for the GL effective charge a(fc 2 ), 

[f^ + ^iL) a(a(//) ' f) = °- (37) 

Solving the equation we may write the GL effective charge as 

^ - (38) 

where C(a) is a scheme independent function. The effective coupling A is defined by 



A(a(u),t) = 7 r-, r , (39) 

-ft (!„« + /*> ^ do -ft)' 

where pi is defined in eq. ( pf) and the integral of 1/ (3(a) is defined in the perturbative form 

a M l 1 /3i 

_-rfo; = — - lnQi(^) + ---. (40) 

/3(a) /waO Po 
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Since 



A{k 2 ) = a{k 2 ) 



l-^-a(k 2 )\na(k 2 ) + --- 

Po 



(41) 



A can be as good an expansion parameter as a. A(k 2 ) is also independent of renormaliza- 
tion scheme. The scheme independence of A(k 2 ) can be easily understood by considering 
the difference of the integral in eq. (|39|) between two schemes. Since the difference is /i 
independent, it cannot depend on a(//) and thus must be a constant. It is in fact given by 



(3'(x) 



Q 0) d x 



7i 



(42) 



where 71 is the first coefficient of the relation between the couplings of the two schemes 



a'(fj,) = a(/i) 1 + 7i«(/i) + 7 2 a(/i) 2 + 



(43) 



The difference in eq. (H2I) is then canceled by the scheme dependence of pi 



Pi =Pi + H, 



(44) 



leading to the scheme independence of A(k 2 ). Note that a similar effective charge was 
considered in relation to the renormalization scheme invariant perturbation by Maximov 
and Vovk [|]. 

Now C(a(fc 2 )) in eq. (EH) may be expanded in an asymptotic series as 



C(a(k 2 )) = d ln(a(A; 2 )) + £ <k [a^ 2 )]^ 1 . 



(45) 



The coefficients Cj can be determined by solving perturbatively a(/c 2 ) in eq. (|3^) with the 
help of eqs. (j39|) , (^0|) in terms of a(/x), and comparing it with the perturbative expansion 



of the GL effective charge in a(fx) using eq. (|T4]). Then it is not difficult to see that q are 
given by the renormalization scheme invariant combinations of the coefficients of the vacuum 
polarization function and (3 function. This expansion turns out to be a critical step for the 



13 



Borel transform of the effective charge. Introduce aS N '(k 2 ) 



a^ N \k 2 ) 



fa + dln(aW(#»)) + £2 r cM N Kk 2 Y- 



(46) 



A{k 2 ) 



Solving a^ N '(k 2 ) in eq. ( ^6|) recursively in ct(yu), it can be seen that this equation generates 
all the higher order renormalon diagrams to order N. 

Let us now consider a modified Borel transform of sS N \k 2 ) defined by 



a ( "'(«M,l) 



x {N \b)db. 



(47) 



Note that this Borel transform of GL effective charge was introduced by Grunberg [|TJ. 
Substituting eq. ([47]) into eq. @, we can write («(//), /i 2 /Q 2 ), which is defined by 
replacing a(A; 2 ) in eq. @ with a^(/c 2 ), as 



cxp 



/3(a) 



{e^ Pl /(6)a^(6)}rf6, (48) 



where 



/(&)= f(t)exp(bf3 \nt)dt. 
Jo 



(49) 



The first IR renormalon singularity arises from the IR divergence in the integral in eq. 
Substituting eq. (|J) into eq. ( fl9"D 



M 



/(&) = / f(t)e b P° lnt dt 



M 



8tt 3 Q 4 Jo 



te bfi \nt dt 



8tt 3 Q 4 2 + b/3 



:i + (2 + 6A))lnM+...), 



(50) 



where M is an arbitrary UV cutoff. Notice that the leading renormalon singularity is cutoff 
independent. 
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The scheme dependence of the Borel transform of D^ N \a(/j,), fj 2 /Q 2 ) is now isolated in 



and 



(51) 



To find the Borel transform explicitly, let us take a renormalization scheme in which the j3 
function is given by the simple form 



/3(a) 



1 — Xa ' 



(52) 



where 



With this (3 function, eq. 
Zhai (BYZ) g, 



Pi 

/V 



(53) 



defines the modified Borel transform by Brown, Yaffe and 



D( N \a,ft 2 /Q-) = I exp 
/o 



f(b)e- hpi e- Xhlnb ^ N \b)db. (54) 



From this we can read off the modified Borel transform of D^ N \a{jj) , /i 2 / Q 2 



D^ ] z (b) = e- bpi f(b)e- Ablnb ^'^(b). 



6 pi f (h\v- Xbhlb *( N )t 



(55) 



Then using the relation between the ordinary Borel transform and the modified one ||, we 
have the ordinary Borel transform 



D W- ltofiQ* 6 6 a ^(l + ^o) 1 " 2 ^' 



(l + 0(2 + 6A))). (56) 



The real problem here is to find the aS N \b). However, with the help of the expansion in 
eq. (H) it is now easy to calculate the Borel transform of a^^A; 2 ). The inverse of eq. ([47]) 
is 

aW(6) = — / e bx eS N \x)dx, (57) 
2m Jc 
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with the contour wrapping around the negative real axis. Here 



1 

x 



A(k 2 ) 
and 



x + ci In aW (x) + E2 Ci[a (N) (x)}^ 1 ' 



Putting 

1 

we can write eq. ([59l) as 

x = y + c 1 Iny-X/^rT: 

i=2 



and eq. fl5T|) as 



where 



2m J - foV 



1 for « = 
Cj = <( Ci for i = 1 
(z — 1) Q for 2 > 2. 

The exponential term in the integrand may be expanded as 

v^iV c 00 ( h\ k k ( N ~ l ) U 



e 



where 



fc=0 ^! ; =fc 1/' 



^ 11, 1 ■ ,+j 



{rii} 
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with the set {n^} of nonnegative integers satisfying 



N-1 N-1 

^2m-i = l, ^2 rii = k. (66) 

i=l 1=1 



Substituting eq. (6|j ) into eq. (^), we finally have 



ubci r oo / h \k k{N-l) 7 TV 



fc=0 "" i=fc ^ i=0 



oo fe(AT-l) TV / -.^7 - 

I 1) llNkltj jk+l+i (c 7 \ 

^ ^ ^ k\T(l + i + l-b Cl ) ' K ' 



k=0 l=k i=0 



This completes the Borel transform of the GL effective charge. 



4 Renormalon residue 

To find the leading renormalon residue of D^ N \b), we have to evaluate b bci a.^ N \b) at the 
first IR renormalon position, b = —2/j3 . If we directly substitute b in eq. (^7|) with b , 
the resulting large order behavior would sum all the contribution from the renormalons to 
order N, but unfortunately this large order behavior does not have a finite limit for iV — ► oo 



TM. The reason for this is that &(b) is singular at the UV and IR renormalon positions, 



and its radius of convergence when it is expanded as in eq. (|67|) is given by the position at 
b — l/A) of the first UV renormalon, which is the closest renormalon to the origin in the 
Borel plane. Therefore we cannot substitute b with bo in eq. ( |67D to correctly evaluate the 
Borel transform at the first IR renormalon. 

This problem can be avoided by introducing an analytic transform of the Borel plane so 
that the closest renormalon to the origin in the new complex plane is the first IR renormalon 
||. Because the singularity of a(6) at the IR renormalon is such that it is finite but has 
divergent derivative [0], we can then express the residue as a convergent series. 

For this purpose, we can take any analytic transform that puts the IR renormalon as the 
closest singularity to the origin, but here we consider a simple form 

z = t=m> (68) 
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with its inverse 

In the 2 — plane, the closest singularity to the origin is the first IR renormalon at 



= -3, (70) 

and all the UV renormalons are pushed beyond z = — 1 on the real axis. It is interesting to 
note that the freedom in choosing this analytic transform is similar to that for the renormal- 
ization scheme. As the renormalization scheme can be optimized for a particular process, the 
analytic transform could be chosen to optimize the perturbative evaluation of the residue. 

Now to find b bci aS N ^(b) at the first IR renormalon, we have to substitute b in eq. (|67D 
with that in eq. ( |69| ) and expand it in Taylor series at z = to order N, and evaluate it at 
z = zq. Thus the Borel transform of GL effective charge at the first IR renormalon is given 
by 



KN = b hc ^ N \b) 



b=b 
N 

E^4 f , (7i) 

M=0 



with 



N fc(JV -i) N N N ( _ l) k +mhNki ^ i 



qM = { M -im-Y, E EEE k \J M ^r } ffr4 +w+J+m ,M, (72) 

fc=o i=k !=oi=om=o K\m\\M - m - ly. 
Here 7]™^ is defined by 

,{«} 



EfV. (73) 



r(n + 1 - x) j=0 

Note that the coefficient is completely determined by sS N \b) in eq. (|67D, and is not 
modified by the Borel transform a^ +m ^(6) of the higher order effective charge. The IR 
residue is then 

00 

«=E^m4'- (74) 

M=0 



18 



It should be emphasized that, though this series is being evaluated at its radius of conver- 
gence, it is convergent because of the fmiteness of a.(b(z)) at z = z . 



5 Large order behavior 

The Borel transform D(b) can now be used in determining the leading large order behavior 
of D(a). Note that the leading large order behavior is determined by the leading Borel 
singularity, and the 1/n correction in the large order behavior corresponds to the 0(b — bo) 
correction in the Borel transform. First, we give the large order behavior in a renormalization 
scheme in which the (3 function is given by eq. (|52|), and then in then sec. 7 discuss a class 
of schemes in which all schemes share a common large order behavior except for a trivial 
scheme dependent term. Using the known result for the vacuum polarization function and 
the j3 function to three-loop and four-loop respectively, we will also give numerical values for 
the large order behavior. 

Let us go back to the Borel transform of D in eq. (|56|). Expanding D^ N \b) at the origin, 
and using the definition of the Borel transform in eq. fl3"S|), the leading large order behavior 
of D( N \a(n),n 2 /Q 2 ) is given by 



e 2 u A 

W = ^_ e -boPi e -Abolnfeo & (Af) / &Q ) n \ n \b 0b -n ^ ^ = -2 (75) 



Then the sequence for the large order behavior mentioned in sec. 1 may be defined as 

K " = Sq^» (76) 

with k n defined in eq (ffH). 

To evaluate the numerical values for k^, we have to find the coefficients q in eq. (^) 
explicitly. The (3 function and the vacuum polarization function in MS scheme to three-loop 
and four-loop respectively are given by JTOj] 
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and 
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where Nf is the number of fermion flavors. Solving a(k 2 ) in eq. 
comparing it with eq. ( |i~4"|) we find 
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(79) 



Note that the use of the vacuum polarization function and the (3 function in MS is allowed, 
because independent of renormalization scheme. 

The kat obtained by substituting eq. (|79|) into eq. (|7l|) is given in Table 1 for several 
flavor numbers. In the table, we see that for a reliable estimation of the large order behavior, 
a higher order calculation beyond the current one is required for the vacuum polarization 
function and the (3 function. 

Not surprisingly, the numbers in the table also suggest that the large Nf limit is the one- 
loop renormalon. This is indeed the case. To see this, note that the following coefficients 
scale as 



0; 



N} 



for i > 0, 
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Table 1: The first three elements of the sequence for the first IR renormalon 
residue in QED. k denotes the residue from the one-loop renor- 
malon. 



q ~ N) 



t-1 



(80) 



for large Nf. Then scaling the variable y in eq. ([62]) by 6 y> it is straightforward to see that 



lim «jy = 1. 



iV f — >oo 



51) 



6 Residue in QCD 

In QCD, there is unfortunately no satisfactory definition of renormalization scheme and 
scale invariant effective charge that may be used in the diagrammatic study of renormalon. 
However, as long as such an effective charge is defined, the formalism developed in QED 
may be used without modification. 

Often in renormalon calculation in QCD, IT(t) in eq. flUD that defines the effective charge 

, for example, as in the 1/Nf approximation combined with 
and the pinch technique |13]|. The pinch technique appears 



is considered in certain limit [p 
"naive-nonabelianization" ||T2| 
to be promising, though presently there is no all-order definition for the effective charge in 
this scheme. In pinch technique, IT(t) at one-loop level is defined by collecting the gluon 
vacuum polarization, and the vacuum polarization-like term in the vertex and box diagrams. 

However, if we are only interested in the calculation of the residue, the definition of the 
effective charge is not required. Indeed the calculation is cunningly simple; it only requires 
the strength of the renormalon singularity and the perturbative calculation of D(a). 

Consider the Borel transform of the current correlation function in QCD. The renormalon 
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singularity of D(b) in QCD 



Di ''> " FW^ (82) 



gives the large order behavior 

a "~7VTTT n!nAb ° 6 o n - (83) 
[A )\ 

To calculate the residue D, consider a function 

R(b) = D{b) (1 -b/b ) 1+Xbo . (84) 
Then because of eq. (|^) , we have 

D = R(b ). (85) 

To avoid the first UV renormalon, we introduce a new variable 2, as we did in QED, which 
is defined by 

1 - /3 b 

with its inverse 

b= fa \l-z 
In the z— plane, the IR renormalon at 



' / 1 ' '87) 



= I (88) 



is the closest singularity to the origin, and so the radius of convergence of the Taylor series 
of D(b(z)) at z = is given by the first IR renormalon. 
Now £) can be expressed in a convergent series form 



b=b 

1+Xbo 



D = D(b) (1 - b/b ) 1+Xbo 
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oo 

= E^o, (89) 
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Table 2: The first three elements of the sequence for the large order behavior 
in QCD. 



where it is straightforward to find r n in terms of the perturbative coefficients a n . Note that 
the series is convergent even if R(b(z)) is not analytic at z = z , because then the radius of 
convergence of the series is given by z = z , and R(b(zo)) is finite. 

Using the perturbative calculation of the current correlation function, and D(a), to three 
loop M, we have 



R(b(z)) 



16 7T 3 



1.333 - 0.748 z - 0.311 z 1 + 0{z A 



(90) 



for Nf = 3. This is in the renormalization scheme in which the one-loop renormalization 



point is same as that of MS scheme, and the {3 function is given in the form in eqs. ([52]), 
([53]). Evaluating this series at the renormalon position at z = z we have 
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(91) 



For several other flavor numbers we give K n in Table 2. 



7 Scheme dependence of large-order behavior 

In section 5, we determined the large order behavior in renormalization schemes in which 
the (3 function is given by the simple form in eq. (0). With this (3 function, the scheme 
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dependence of the large order behavior arises only through the factor 



e- baPl (92) 

in eq. (|76|) . In fact, this result is more general. All renormalization schemes for which the 
coefficients of the (3 function do not grow faster than a n share the same large-large order 
behavior except for the scheme dependence in eq. (|92|). 

To see this, let us consider two renormalization schemes (say, unprimed and primed ) in 
which the relation between the couplings is given by eq. (pE3|). We now assume that the large 
order behavior by the first UV renormalon is extracted out so that the leading large order 
behavior of a scheme independent subamplitude D' of D is given by the first IR renormalon. 
Then using 

oo oo 

D'=J2 a„«(/i)" +1 = £ «;«'(/i) n+1 , (93) 

n=0 n=0 

we get the relation between a n and a' n 

a „ = <(l + (n+ l)%l 7l + !^±i)%l 7 f + + (94) 

V a n 1 a n a n J 

With the large order behavior of a' n in the form in eq. flU), this equation becomes 

a' n = a n e- b0 ' n (1 + 0(1 /n)) (95) 

provided 

lim = const. (96) 

Thus if 7„ does not grow faster than 

(n-l)!n% n , (97) 



the scheme dependence of large order behavior is given by the simple relation in eq. flOg) . 
In fact the large order behavior we found in sec. 5 exactly transforms according to eq. ( |95"D 
under scheme changes. 

We may now translate the limit in eq. (pB) on 7„ to that of (3 n . Let the (3 function in 
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the unprimed scheme is given by the simple form in eq. (52) . Then eq. (|43|) gives the f3'(a!) 
in the form 



where 



n=0 n=0 

p a 2 



J2(n + l) ln a n 
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n=0 



0'n = E^" fc (^ + l)7 fc 



fc=0 



~ 7i7 n (l + 0(l/n)) (99) 
Inverting eq. ( fSj) to express a in terms of a', and substituting it in eq. ([98D , we have 

= /3' n (l + 0(l/n)) = n 7 „(l + 0(l/n)). (100) 



Then the restriction on 7n in eq. (|97D implies that any renormalization scheme in which /3 n 
does not grow faster than 

n\n u b n (101) 

has the same large order behavior (except for the factor eq. (j92[) ) as in a scheme for which 
the (3 function is given by eq. (|52|) . 
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